3.1.1 \(\int x^4 \text {sech}^{-1}(a x)^2 \, dx\) [1]

3.1.1.1 Optimal result
3.1.1.2 Mathematica [A] (verified)
3.1.1.3 Rubi [A] (verified)
3.1.1.4 Maple [A] (verified)
3.1.1.5 Fricas [F]
3.1.1.6 Sympy [F]
3.1.1.7 Maxima [F(-1)]
3.1.1.8 Giac [F]
3.1.1.9 Mupad [F(-1)]

3.1.1.1 Optimal result

Integrand size = 10, antiderivative size = 164 \[ \int x^4 \text {sech}^{-1}(a x)^2 \, dx=-\frac {3 x}{20 a^4}-\frac {x^3}{30 a^2}-\frac {3 x \sqrt {\frac {1-a x}{1+a x}} (1+a x) \text {sech}^{-1}(a x)}{20 a^4}-\frac {x^3 \sqrt {\frac {1-a x}{1+a x}} (1+a x) \text {sech}^{-1}(a x)}{10 a^2}+\frac {1}{5} x^5 \text {sech}^{-1}(a x)^2-\frac {3 \text {sech}^{-1}(a x) \arctan \left (e^{\text {sech}^{-1}(a x)}\right )}{10 a^5}+\frac {3 i \operatorname {PolyLog}\left (2,-i e^{\text {sech}^{-1}(a x)}\right )}{20 a^5}-\frac {3 i \operatorname {PolyLog}\left (2,i e^{\text {sech}^{-1}(a x)}\right )}{20 a^5} \]

output
-3/20*x/a^4-1/30*x^3/a^2+1/5*x^5*arcsech(a*x)^2-3/10*arcsech(a*x)*arctan(1 
/a/x+(1/a/x-1)^(1/2)*(1+1/a/x)^(1/2))/a^5+3/20*I*polylog(2,-I*(1/a/x+(1/a/ 
x-1)^(1/2)*(1+1/a/x)^(1/2)))/a^5-3/20*I*polylog(2,I*(1/a/x+(1/a/x-1)^(1/2) 
*(1+1/a/x)^(1/2)))/a^5-3/20*x*(a*x+1)*arcsech(a*x)*((-a*x+1)/(a*x+1))^(1/2 
)/a^4-1/10*x^3*(a*x+1)*arcsech(a*x)*((-a*x+1)/(a*x+1))^(1/2)/a^2
 
3.1.1.2 Mathematica [A] (verified)

Time = 0.42 (sec) , antiderivative size = 182, normalized size of antiderivative = 1.11 \[ \int x^4 \text {sech}^{-1}(a x)^2 \, dx=\frac {-9 a x-2 a^3 x^3-9 a x \sqrt {\frac {1-a x}{1+a x}} (1+a x) \text {sech}^{-1}(a x)-6 a^3 x^3 \sqrt {\frac {1-a x}{1+a x}} (1+a x) \text {sech}^{-1}(a x)+12 a^5 x^5 \text {sech}^{-1}(a x)^2+9 i \text {sech}^{-1}(a x) \log \left (1-i e^{-\text {sech}^{-1}(a x)}\right )-9 i \text {sech}^{-1}(a x) \log \left (1+i e^{-\text {sech}^{-1}(a x)}\right )+9 i \operatorname {PolyLog}\left (2,-i e^{-\text {sech}^{-1}(a x)}\right )-9 i \operatorname {PolyLog}\left (2,i e^{-\text {sech}^{-1}(a x)}\right )}{60 a^5} \]

input
Integrate[x^4*ArcSech[a*x]^2,x]
 
output
(-9*a*x - 2*a^3*x^3 - 9*a*x*Sqrt[(1 - a*x)/(1 + a*x)]*(1 + a*x)*ArcSech[a* 
x] - 6*a^3*x^3*Sqrt[(1 - a*x)/(1 + a*x)]*(1 + a*x)*ArcSech[a*x] + 12*a^5*x 
^5*ArcSech[a*x]^2 + (9*I)*ArcSech[a*x]*Log[1 - I/E^ArcSech[a*x]] - (9*I)*A 
rcSech[a*x]*Log[1 + I/E^ArcSech[a*x]] + (9*I)*PolyLog[2, (-I)/E^ArcSech[a* 
x]] - (9*I)*PolyLog[2, I/E^ArcSech[a*x]])/(60*a^5)
 
3.1.1.3 Rubi [A] (verified)

Time = 0.59 (sec) , antiderivative size = 168, normalized size of antiderivative = 1.02, number of steps used = 11, number of rules used = 10, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 1.000, Rules used = {6839, 5941, 3042, 4673, 3042, 4673, 3042, 4668, 2715, 2838}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int x^4 \text {sech}^{-1}(a x)^2 \, dx\)

\(\Big \downarrow \) 6839

\(\displaystyle -\frac {\int a^5 x^5 \sqrt {\frac {1-a x}{a x+1}} (a x+1) \text {sech}^{-1}(a x)^2d\text {sech}^{-1}(a x)}{a^5}\)

\(\Big \downarrow \) 5941

\(\displaystyle -\frac {\frac {2}{5} \int a^5 x^5 \text {sech}^{-1}(a x)d\text {sech}^{-1}(a x)-\frac {1}{5} a^5 x^5 \text {sech}^{-1}(a x)^2}{a^5}\)

\(\Big \downarrow \) 3042

\(\displaystyle -\frac {-\frac {1}{5} a^5 x^5 \text {sech}^{-1}(a x)^2+\frac {2}{5} \int \text {sech}^{-1}(a x) \csc \left (i \text {sech}^{-1}(a x)+\frac {\pi }{2}\right )^5d\text {sech}^{-1}(a x)}{a^5}\)

\(\Big \downarrow \) 4673

\(\displaystyle -\frac {\frac {2}{5} \left (\frac {3}{4} \int a^3 x^3 \text {sech}^{-1}(a x)d\text {sech}^{-1}(a x)+\frac {a^3 x^3}{12}+\frac {1}{4} a^3 x^3 \sqrt {\frac {1-a x}{a x+1}} (a x+1) \text {sech}^{-1}(a x)\right )-\frac {1}{5} a^5 x^5 \text {sech}^{-1}(a x)^2}{a^5}\)

\(\Big \downarrow \) 3042

\(\displaystyle -\frac {-\frac {1}{5} a^5 x^5 \text {sech}^{-1}(a x)^2+\frac {2}{5} \left (\frac {3}{4} \int \text {sech}^{-1}(a x) \csc \left (i \text {sech}^{-1}(a x)+\frac {\pi }{2}\right )^3d\text {sech}^{-1}(a x)+\frac {a^3 x^3}{12}+\frac {1}{4} a^3 x^3 \sqrt {\frac {1-a x}{a x+1}} (a x+1) \text {sech}^{-1}(a x)\right )}{a^5}\)

\(\Big \downarrow \) 4673

\(\displaystyle -\frac {\frac {2}{5} \left (\frac {3}{4} \left (\frac {1}{2} \int a x \text {sech}^{-1}(a x)d\text {sech}^{-1}(a x)+\frac {a x}{2}+\frac {1}{2} a x \sqrt {\frac {1-a x}{a x+1}} (a x+1) \text {sech}^{-1}(a x)\right )+\frac {a^3 x^3}{12}+\frac {1}{4} a^3 x^3 \sqrt {\frac {1-a x}{a x+1}} (a x+1) \text {sech}^{-1}(a x)\right )-\frac {1}{5} a^5 x^5 \text {sech}^{-1}(a x)^2}{a^5}\)

\(\Big \downarrow \) 3042

\(\displaystyle -\frac {-\frac {1}{5} a^5 x^5 \text {sech}^{-1}(a x)^2+\frac {2}{5} \left (\frac {3}{4} \left (\frac {1}{2} \int \text {sech}^{-1}(a x) \csc \left (i \text {sech}^{-1}(a x)+\frac {\pi }{2}\right )d\text {sech}^{-1}(a x)+\frac {a x}{2}+\frac {1}{2} a x \sqrt {\frac {1-a x}{a x+1}} (a x+1) \text {sech}^{-1}(a x)\right )+\frac {a^3 x^3}{12}+\frac {1}{4} a^3 x^3 \sqrt {\frac {1-a x}{a x+1}} (a x+1) \text {sech}^{-1}(a x)\right )}{a^5}\)

\(\Big \downarrow \) 4668

\(\displaystyle -\frac {-\frac {1}{5} a^5 x^5 \text {sech}^{-1}(a x)^2+\frac {2}{5} \left (\frac {3}{4} \left (\frac {1}{2} \left (-i \int \log \left (1-i e^{\text {sech}^{-1}(a x)}\right )d\text {sech}^{-1}(a x)+i \int \log \left (1+i e^{\text {sech}^{-1}(a x)}\right )d\text {sech}^{-1}(a x)+2 \text {sech}^{-1}(a x) \arctan \left (e^{\text {sech}^{-1}(a x)}\right )\right )+\frac {a x}{2}+\frac {1}{2} a x \sqrt {\frac {1-a x}{a x+1}} (a x+1) \text {sech}^{-1}(a x)\right )+\frac {a^3 x^3}{12}+\frac {1}{4} a^3 x^3 \sqrt {\frac {1-a x}{a x+1}} (a x+1) \text {sech}^{-1}(a x)\right )}{a^5}\)

\(\Big \downarrow \) 2715

\(\displaystyle -\frac {-\frac {1}{5} a^5 x^5 \text {sech}^{-1}(a x)^2+\frac {2}{5} \left (\frac {3}{4} \left (\frac {1}{2} \left (-i \int e^{-\text {sech}^{-1}(a x)} \log \left (1-i e^{\text {sech}^{-1}(a x)}\right )de^{\text {sech}^{-1}(a x)}+i \int e^{-\text {sech}^{-1}(a x)} \log \left (1+i e^{\text {sech}^{-1}(a x)}\right )de^{\text {sech}^{-1}(a x)}+2 \text {sech}^{-1}(a x) \arctan \left (e^{\text {sech}^{-1}(a x)}\right )\right )+\frac {a x}{2}+\frac {1}{2} a x \sqrt {\frac {1-a x}{a x+1}} (a x+1) \text {sech}^{-1}(a x)\right )+\frac {a^3 x^3}{12}+\frac {1}{4} a^3 x^3 \sqrt {\frac {1-a x}{a x+1}} (a x+1) \text {sech}^{-1}(a x)\right )}{a^5}\)

\(\Big \downarrow \) 2838

\(\displaystyle -\frac {-\frac {1}{5} a^5 x^5 \text {sech}^{-1}(a x)^2+\frac {2}{5} \left (\frac {a^3 x^3}{12}+\frac {1}{4} a^3 x^3 \sqrt {\frac {1-a x}{a x+1}} (a x+1) \text {sech}^{-1}(a x)+\frac {3}{4} \left (\frac {1}{2} \left (2 \text {sech}^{-1}(a x) \arctan \left (e^{\text {sech}^{-1}(a x)}\right )-i \operatorname {PolyLog}\left (2,-i e^{\text {sech}^{-1}(a x)}\right )+i \operatorname {PolyLog}\left (2,i e^{\text {sech}^{-1}(a x)}\right )\right )+\frac {a x}{2}+\frac {1}{2} a x \sqrt {\frac {1-a x}{a x+1}} (a x+1) \text {sech}^{-1}(a x)\right )\right )}{a^5}\)

input
Int[x^4*ArcSech[a*x]^2,x]
 
output
-((-1/5*(a^5*x^5*ArcSech[a*x]^2) + (2*((a^3*x^3)/12 + (a^3*x^3*Sqrt[(1 - a 
*x)/(1 + a*x)]*(1 + a*x)*ArcSech[a*x])/4 + (3*((a*x)/2 + (a*x*Sqrt[(1 - a* 
x)/(1 + a*x)]*(1 + a*x)*ArcSech[a*x])/2 + (2*ArcSech[a*x]*ArcTan[E^ArcSech 
[a*x]] - I*PolyLog[2, (-I)*E^ArcSech[a*x]] + I*PolyLog[2, I*E^ArcSech[a*x] 
])/2))/4))/5)/a^5)
 

3.1.1.3.1 Defintions of rubi rules used

rule 2715
Int[Log[(a_) + (b_.)*((F_)^((e_.)*((c_.) + (d_.)*(x_))))^(n_.)], x_Symbol] 
:> Simp[1/(d*e*n*Log[F])   Subst[Int[Log[a + b*x]/x, x], x, (F^(e*(c + d*x) 
))^n], x] /; FreeQ[{F, a, b, c, d, e, n}, x] && GtQ[a, 0]
 

rule 2838
Int[Log[(c_.)*((d_) + (e_.)*(x_)^(n_.))]/(x_), x_Symbol] :> Simp[-PolyLog[2 
, (-c)*e*x^n]/n, x] /; FreeQ[{c, d, e, n}, x] && EqQ[c*d, 1]
 

rule 3042
Int[u_, x_Symbol] :> Int[DeactivateTrig[u, x], x] /; FunctionOfTrigOfLinear 
Q[u, x]
 

rule 4668
Int[csc[(e_.) + Pi*(k_.) + (Complex[0, fz_])*(f_.)*(x_)]*((c_.) + (d_.)*(x_ 
))^(m_.), x_Symbol] :> Simp[-2*(c + d*x)^m*(ArcTanh[E^((-I)*e + f*fz*x)/E^( 
I*k*Pi)]/(f*fz*I)), x] + (-Simp[d*(m/(f*fz*I))   Int[(c + d*x)^(m - 1)*Log[ 
1 - E^((-I)*e + f*fz*x)/E^(I*k*Pi)], x], x] + Simp[d*(m/(f*fz*I))   Int[(c 
+ d*x)^(m - 1)*Log[1 + E^((-I)*e + f*fz*x)/E^(I*k*Pi)], x], x]) /; FreeQ[{c 
, d, e, f, fz}, x] && IntegerQ[2*k] && IGtQ[m, 0]
 

rule 4673
Int[(csc[(e_.) + (f_.)*(x_)]*(b_.))^(n_)*((c_.) + (d_.)*(x_)), x_Symbol] :> 
 Simp[(-b^2)*(c + d*x)*Cot[e + f*x]*((b*Csc[e + f*x])^(n - 2)/(f*(n - 1))), 
 x] + (-Simp[b^2*d*((b*Csc[e + f*x])^(n - 2)/(f^2*(n - 1)*(n - 2))), x] + S 
imp[b^2*((n - 2)/(n - 1))   Int[(c + d*x)*(b*Csc[e + f*x])^(n - 2), x], x]) 
 /; FreeQ[{b, c, d, e, f}, x] && GtQ[n, 1] && NeQ[n, 2]
 

rule 5941
Int[(x_)^(m_.)*Sech[(a_.) + (b_.)*(x_)^(n_.)]^(p_.)*Tanh[(a_.) + (b_.)*(x_) 
^(n_.)]^(q_.), x_Symbol] :> Simp[(-x^(m - n + 1))*(Sech[a + b*x^n]^p/(b*n*p 
)), x] + Simp[(m - n + 1)/(b*n*p)   Int[x^(m - n)*Sech[a + b*x^n]^p, x], x] 
 /; FreeQ[{a, b, p}, x] && RationalQ[m] && IntegerQ[n] && GeQ[m - n, 0] && 
EqQ[q, 1]
 

rule 6839
Int[((a_.) + ArcSech[(c_.)*(x_)]*(b_.))^(n_)*(x_)^(m_.), x_Symbol] :> Simp[ 
-(c^(m + 1))^(-1)   Subst[Int[(a + b*x)^n*Sech[x]^(m + 1)*Tanh[x], x], x, A 
rcSech[c*x]], x] /; FreeQ[{a, b, c}, x] && IntegerQ[n] && IntegerQ[m] && (G 
tQ[n, 0] || LtQ[m, -1])
 
3.1.1.4 Maple [A] (verified)

Time = 0.62 (sec) , antiderivative size = 280, normalized size of antiderivative = 1.71

method result size
derivativedivides \(\frac {\frac {\left (-6 \,\operatorname {arcsech}\left (a x \right ) \sqrt {-\frac {a x -1}{a x}}\, \sqrt {\frac {a x +1}{a x}}\, a^{3} x^{3}+12 a^{4} x^{4} \operatorname {arcsech}\left (a x \right )^{2}-9 \,\operatorname {arcsech}\left (a x \right ) \sqrt {-\frac {a x -1}{a x}}\, \sqrt {\frac {a x +1}{a x}}\, a x -2 a^{2} x^{2}-9\right ) a x}{60}+\frac {3 i \operatorname {arcsech}\left (a x \right ) \ln \left (1+i \left (\frac {1}{a x}+\sqrt {\frac {1}{a x}-1}\, \sqrt {1+\frac {1}{a x}}\right )\right )}{20}-\frac {3 i \operatorname {arcsech}\left (a x \right ) \ln \left (1-i \left (\frac {1}{a x}+\sqrt {\frac {1}{a x}-1}\, \sqrt {1+\frac {1}{a x}}\right )\right )}{20}+\frac {3 i \operatorname {dilog}\left (1+i \left (\frac {1}{a x}+\sqrt {\frac {1}{a x}-1}\, \sqrt {1+\frac {1}{a x}}\right )\right )}{20}-\frac {3 i \operatorname {dilog}\left (1-i \left (\frac {1}{a x}+\sqrt {\frac {1}{a x}-1}\, \sqrt {1+\frac {1}{a x}}\right )\right )}{20}}{a^{5}}\) \(280\)
default \(\frac {\frac {\left (-6 \,\operatorname {arcsech}\left (a x \right ) \sqrt {-\frac {a x -1}{a x}}\, \sqrt {\frac {a x +1}{a x}}\, a^{3} x^{3}+12 a^{4} x^{4} \operatorname {arcsech}\left (a x \right )^{2}-9 \,\operatorname {arcsech}\left (a x \right ) \sqrt {-\frac {a x -1}{a x}}\, \sqrt {\frac {a x +1}{a x}}\, a x -2 a^{2} x^{2}-9\right ) a x}{60}+\frac {3 i \operatorname {arcsech}\left (a x \right ) \ln \left (1+i \left (\frac {1}{a x}+\sqrt {\frac {1}{a x}-1}\, \sqrt {1+\frac {1}{a x}}\right )\right )}{20}-\frac {3 i \operatorname {arcsech}\left (a x \right ) \ln \left (1-i \left (\frac {1}{a x}+\sqrt {\frac {1}{a x}-1}\, \sqrt {1+\frac {1}{a x}}\right )\right )}{20}+\frac {3 i \operatorname {dilog}\left (1+i \left (\frac {1}{a x}+\sqrt {\frac {1}{a x}-1}\, \sqrt {1+\frac {1}{a x}}\right )\right )}{20}-\frac {3 i \operatorname {dilog}\left (1-i \left (\frac {1}{a x}+\sqrt {\frac {1}{a x}-1}\, \sqrt {1+\frac {1}{a x}}\right )\right )}{20}}{a^{5}}\) \(280\)

input
int(x^4*arcsech(a*x)^2,x,method=_RETURNVERBOSE)
 
output
1/a^5*(1/60*(-6*arcsech(a*x)*(-(a*x-1)/a/x)^(1/2)*((a*x+1)/a/x)^(1/2)*a^3* 
x^3+12*a^4*x^4*arcsech(a*x)^2-9*arcsech(a*x)*(-(a*x-1)/a/x)^(1/2)*((a*x+1) 
/a/x)^(1/2)*a*x-2*a^2*x^2-9)*a*x+3/20*I*arcsech(a*x)*ln(1+I*(1/a/x+(1/a/x- 
1)^(1/2)*(1+1/a/x)^(1/2)))-3/20*I*arcsech(a*x)*ln(1-I*(1/a/x+(1/a/x-1)^(1/ 
2)*(1+1/a/x)^(1/2)))+3/20*I*dilog(1+I*(1/a/x+(1/a/x-1)^(1/2)*(1+1/a/x)^(1/ 
2)))-3/20*I*dilog(1-I*(1/a/x+(1/a/x-1)^(1/2)*(1+1/a/x)^(1/2))))
 
3.1.1.5 Fricas [F]

\[ \int x^4 \text {sech}^{-1}(a x)^2 \, dx=\int { x^{4} \operatorname {arsech}\left (a x\right )^{2} \,d x } \]

input
integrate(x^4*arcsech(a*x)^2,x, algorithm="fricas")
 
output
integral(x^4*arcsech(a*x)^2, x)
 
3.1.1.6 Sympy [F]

\[ \int x^4 \text {sech}^{-1}(a x)^2 \, dx=\int x^{4} \operatorname {asech}^{2}{\left (a x \right )}\, dx \]

input
integrate(x**4*asech(a*x)**2,x)
 
output
Integral(x**4*asech(a*x)**2, x)
 
3.1.1.7 Maxima [F(-1)]

Timed out. \[ \int x^4 \text {sech}^{-1}(a x)^2 \, dx=\text {Timed out} \]

input
integrate(x^4*arcsech(a*x)^2,x, algorithm="maxima")
 
output
Timed out
 
3.1.1.8 Giac [F]

\[ \int x^4 \text {sech}^{-1}(a x)^2 \, dx=\int { x^{4} \operatorname {arsech}\left (a x\right )^{2} \,d x } \]

input
integrate(x^4*arcsech(a*x)^2,x, algorithm="giac")
 
output
integrate(x^4*arcsech(a*x)^2, x)
 
3.1.1.9 Mupad [F(-1)]

Timed out. \[ \int x^4 \text {sech}^{-1}(a x)^2 \, dx=\int x^4\,{\mathrm {acosh}\left (\frac {1}{a\,x}\right )}^2 \,d x \]

input
int(x^4*acosh(1/(a*x))^2,x)
 
output
int(x^4*acosh(1/(a*x))^2, x)